Density Profiles of Strongly Interacting Trapped Fermi Gases 
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We study density profiles in trapped fermionic gases, near Feshbach resonances, at all T < Tc and in the 
near-BEC and unitary regimes. For the latter, we characterize and quantify the generally neglected contribution 
from noncondensed Cooper pairs. As a consequence of these pairs, our profiles are rather well fit to a Thomas- 
Fermi (TF) functional form, and equally well fit to experimental data. Our work lends support to the notion that 
TF fits can be used in an experimental context to obtain information about the temperature. 



There is now strong, but not unambiguous, evidence that 
fermionic superfluidity has been observed 

Gliiilli 

in trapped atomic gases. What is particularly exciting about 
these systems is that the strength of the pairing interaction 
can be arbitrarily tuned (via Feshbach resonance effects) so 
that superfluidity occurs in a regime which goes beyond tradi- 
tional weak coupling BCS theory, towards Bose-Einstein con- 
densation (BEC) of pre-formed pairs. These experiments are 
revealing the nature of superconductivity and superfluidity in 
a hitherto unexplored regime. 

Experimental proof of superfluidity is not straightforward, 
except in the BEC regime where the particle density profiles 
acquire a bi-modal form. This bi-modality is similar to obser- 
vations in bosonic atomic systems, which have been system- 
atically studied for the past decade tZJ. In the intermediate 
regime between BCS and BEC, the measured density profiles 
1 8, "q"] in traps do not contain any obvious signatures of the 
condensate. Nevertheless, there is a significant body of cir- 
cumstantial evidence in support of superfluidity in this regime. 
This support is based on fast sweep experiments Ol|^, col- 
lective mode measurements \2s A, as well as pairing gap 
measurements using radio frequency (RF) techniques O. 

The goals of this paper are to compute the particle density 
profiles at all T < Tc (where Tc is the superfluid transition 
temperature) in the near-BEC and unitary regimes and to ad- 
dress the implications for experiment. For the former our cal- 
culation of the contribution from noncondensed bosons leads 
to important corrections to estimates of the condensate frac- 
tion based on the conventional Gaussian form. For the uni- 
tary case, these non-condensed pairs smooth out the, other- 
wise abrupt, transition between the condensate and the ther- 
mal background of fermionic excitations; this explains why 
the measured density profiles appear to be so featureless . 

The unitary or strongly interacting Fermi gas has been the 
focus of attention by the community. Some time ago it was 
found fi*] that the profiles were reasonably well described by 
a Thomas-Fermi (TF) functional form at zero T, and in re- 
cent work this procedure has been extended to finite tem- 
peratures. Importantly, there has been no particular theoret- 
ical support for these TF fits. Below Tc all previous theo- 
retical work (which either ignored non-condensed pair states 
flli [l^ or used a different ground state fl^) has predicted 
strong deviations from this T-dependent TF functional form. 
It is extremely important, thus, to have a better understanding 
of the spatial profiles. Our work lends theoretical support to 



the viewpoint fioll that TF fits (with proper calibration) can be 
used in an experimental context to obtain information about 
the temperature. In addition to establishing these TF fits to 
theory, we compare theory and experiment directly via one 
dimensional representations of their respective profiles. We 
demonstrate remarkable agreement in the shapes of the two 
profiles. One can infer from this comparison, experimental 
support, albeit indirect, for the presence of a superfluid con- 
densate. 

The first generation theories 1 14] have focused on a BCS- 
like ground state (l^ which is readily generalized to ac- 
commodate a 'two channel" variant y^] in which there are 
Feshbach bosonic (FB) degrees of freedom present as well. 
This approach has met with some initial success in addressing 
collective mode experiments fl7i ITSi] and pairing-gap spec- 
troscopy fl^ll . The excitations of this standard ground state 
] 14, 16] consist of two types: noncondensed fermion pairs 
hybridized with FB, and Bogoliubov-like fermionic excita- 
tions having = (et — /i)^ + A2(T), where et is the 
free fermion kinetic energy. Note that the "gap parameter" 
A(T) is in general different from the superconducting order 
parameter, Asc- Recent RF experiments 0] have been an- 
alyzed \Vf\ using our formalism to suggest that the Bogoli- 
ubov quasi-particles have an energy gap or pseudogap, which 
is present well above Tc, and, thus, not directly related to the 
order parameter Additional support for this "pseudogap" has 
recently been reported in another very different class of exper- 
iments ]20]. 

We summarize the self-consistent equations iilll . in the 
presence of a spherical trap, treated at the level of the lo- 
cal density approximation (LDA) with trap potential V(r) = 
^mw^r^. Tc is defined as the highest temperature at which 
the self-consistent equations are satisfied precisely at the cen- 
ter At a temperature T lower than Tc the superfluid region 
extends to a finite radius Rac- The particles outside this ra- 
dius are in a normal state, with or without a pseudogap. Our 
self-consistent equations are given in terms of the Feshbach 
coupling constant g and inter-fermion attractive interaction U 
by a gap equation 



E 



1 - 2/(^k) 
2£;w 







(1) 



which is imposed only when yu^air ir) = 0. The pseudogap 
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contribution to A^{T) = A^^iT) + A^giT) is given by 
q 

The various residues, Z and Z},, which are of no particular 
physical importance here, are described in Ref.il4j The quan- 
tity ilq is the pair dispersion, and fipair and Hboson are the 
effective chemical potentials of the pairs and FB. 

We introduce units fi = \, fermion mass m — ^, Fermi 
momentum kp = 1, and Fermi energy 

In the figures below, we express length in units of the Thomas- 
Fermi radius Rtf = y/2EF/{muj^) = 2i3N)^/^ /kp; the 
density n{r) and total particle number N = J d^rn{r) are 
normalized by kp and {kFRrp)^, respectively. The density 
of particles at radius r can be written as 

n(r) = 2n|J + — ^ 6(rig - fiboson) 

q 

+ 2j2H{l-f{E^))+<f{E,,)] , (3) 

k 

where = g'^Al^/[{iy -2fi + 2V{r))U + g'^f is the molec- 
ular Bose condensate which is only present for r < Rsc- The 
important chemical potential iipair = Hboson i22il is identi- 
cally zero in the superfluid region r < Rsc , and must be 
solved for self-consistently at larger radii. Our calculations 
proceed by solving numerically the self-consistent equations. 

We now decompose the density profiles into components 
associated with the condensate, the noncondensed pairs and 
the fermions. This raises a central issue of this paper What 
is the most suitable definition of the condensate ratio in a 
fermionic superfluid? Indeed, in strict BCS theory there are 
two alternatives |23], one based on the superfluid density, and 
the other based on associating the condensate with the per- 
turbation of the superfluid state relative to the underlying /ree 
Fermi gas. In BCS theory with this second definition, the zero 
temperature condensate fraction is of order Td Ef, far from 
the value 100%, which one obtains from the superfluid den- 
sity. Physically, this second definition reflects the fraction of 
the original Fermi liquid states which are modified substan- 
tially by pairing. We explore both of these here, since they 
are expected to enter in different physical contexts. In the 
BEC regime (where fermions are absent) there is no distinc- 
tion between the two decompositions, and thus the condensate 
fraction is uniquely defined. 

In the present approach to the BCS-BEC crossover picture 
it is relatively straightforward to deduce fT?] the (local) super- 
fluid density, = rts(r). For the one-channel model we find 
lll4ll a simple result for Us, as well as for the fermionic quasi- 
particle (ngp) and pair contributions {Upair) to the difference 
n - Us-. 

= ^nf^^(A), (4a) 
npa.r = ^nf ^^(A) = nf ^^(A) - n, , (4b) 
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FIG. 1 : (color online) Density profile and its decomposition in the 
near-BEC limit (l/kpa = 3.0) at T/Tc = 0.5. Shown are the actual 
plots on the left for the total n(r), the noncondensed pairs Upair, and 
the condensate Us . The right panel shows how these two (indicated 
by the shaded areas from top to bottom) add together to give n{r). 

UQp = n-nf^^iA), (4c) 

where nf^'^ is the usual superfluid density for a gas of 
fermions within BCS theory, except that here the full excita- 
tion gap A appears instead of the order parameter. In the one- 
channel case the order parameter Age is equal to the Cooper 
condensate contribution (which we call Age)- For the broad 
Feshbach resonances of ^''K and ^Li, these one-channel re- 
sults are appropriate for the unitary and BCS regimes of the 
two-channel case. In the BEC regime (/i < 0), the con- 
tribution of the fermionic quasiparticles becomes negligible. 
Therefore, the one-channel equations will be reduced to 

AL ^pg 

Interestingly, this result is also valid for the two-channel prob- 
lem in the BEC regime. 

Formally, we can address the second decomposition of 
the density profiles by writing the particle density in terms 
of the single particle Green's function G{K) = Go{K) + 
G{){K)Yj{K)G{K) where Go represents free fermions, and 
K = {iujn , k) is the four-momentum. The second term is 
then further split into the condensed (oc A^^) and noncon- 
densed (oc Apg) components. Summing over the Matsubara 
frequencies ojn, and adding the FB contribution one obtains 

fis = 2ZAl^, (6a) 

flpair = 2 ^ b{nq - fipair) 2ZAp^ , (6b) 

f^flZ = 2^/(6u-/i + VM), (6c) 

k 

which correspond to condensate density, finite momentum 
pair/boson density, and free fermion density, respectively. 

We characterize the crossover regime by the parameter 
1/kpa where a is the inter-atomic s-wave scattering length 
and kp, the Fermi momentum. A typical density profile for 
the near-BEC limit with 1/kpa = 3.0 is shown in Fig. [0 
for T/Tc = 0.5. This plot is in the physically accessible 
near-BEC regime where the bosons are still primarily Cooper 
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FIG. 2: (color online) Density profiles and their different decompo- 
sitions for the unitary case (l/kpa — 0) atT — 0.75Tc. The upper 
and lower panels correspond to Eqs. j6| and Eqs. J4}, respectively. 
On the left are shown the actual plots for the total, the noncondensed 
pairs, condensate, and fermion density. Shown on the right is how 
the last three, as indicated by the shaded areas from top to bottom of 
each panel, add up to the local total density n{r). 



pairs, as distinguished from Feshbach bosons. Based on ex- 
periments f3, 5] and related theory this is not yet 
in the deep BEC where the collective modes are expected to 
exhibit the characteristics of true bosons jjjt . Thus one ex- 
pects the ground state wavefunction to be applicable here. The 
left panel corresponds to actual plots of the condensed (green 
curve) and noncondensed pair contributions (red) while the 
shaded regions on the right represent from top to bottom the 
condensed (green) and noncondensed bosons (red), respec- 
tively. For this BEC plot, there are no fermions. Importantly, 
the contribution from the noncondensed bosons is essentially 
constant in r until Rsc, reflecting the fact that they have a gap- 
less excitation spectrum (fipair = 0). However, once /ipair is 
non-zero at the trap edge the number of noncondensed bosons 
(upatr) drops rapidly. 

Fig. [2is consistent with the experimentally observed pro- 
file shapes but it can be contrasted with other theoretical 
results in the literature which predict non-monotonic features 
llisll . Because all direct interactions are via fermions, the char- 
acter of the BEC profile is different from that of true interact- 
ing bosons. These differences also appear in the context of 
collective mode experiments lll7lll8ll25ll . As for true bosons 
f?,'2^, the constraint that the bosons are gapless in the super- 
fluid region is important for determining their density distribu- 
tion. A clear bi-modal feature or "kink" at Rsc is present and 
provides evidence for the existence of a condensate. It should 
be noted that the inferred fraction of the condensate may be 
significantly smaller than found here - if a Gaussian form is 
assumed throughout the trap (as is the experimental conven- 
tion) for estimating the contribution from noncondensed pairs. 

In the upper and lower panels of Fig.|2lwe plot the density 
distributions in the unitary limit, 1/kpa = 0, as a function 
of radius r at T/Tc = 0.75, and show their different de- 



compositions following Eqs. (|6} and Eqs. (0), respectively. 
The superfluid-based decomposition (lower panels) appears 
to be more relevant to thermodynamics and recent RF exper- 
iments ^ since it incorporates the excitation gap of the 
fermions. The decomposition in the upper panels, based on 
free fermions, was studied in detail in Ref. 13, for a somewhat 
different ground state. 

The shaded regions on the right of each figure correspond 
(from top to bottom) to density profiles of noncondensed 
pairs, condensate and fermionic quasiparticles. In both de- 
compositions the noncondensed pairs play a significant role. 
It can be seen from the plots on the left of Fig.|2]that the non- 
condensed pairs show a relatively flat density distribution for 
r < Rsc, just as in the near-BEC limit of Fig. [0 This be- 
havior similarly arises from the vanishing of ^pair [r) in the 
superfluid region. 

One can see that the condensate has substantially more 
weight for the n^-based decomposition [Eq. 0, Fig.|2] lower 
panels]. Conversely the excited fermions for this case are sig- 
nificantly less important, since they experience a large gap in 
their excitation spectrum. The contribution of these fermions 
is concentrated at the trap edge. Indeed, for this case, because 
they represent two sides of the same coin, the fermions and 
noncondensed bosons generally appear in different regions of 
the trap. The latter are associated with the regions where the 
excitation gap is largest (and the fermion states are quasi- 
bound into "bosonic" like states), and the former are found 
where the gap is smallest. 

In Fig.|3]we compare theory (at T/Tp — 0.19) and experi- 
ment for the unitary case. The profiles shown are well within 
the superfluid phase {Tc ~ 0.3Tf at unitarity). This figure 
presents Thomas-Fermi fits lHoll to the experimental ((3^) and 
theoretical (|3j)) profiles as well as their comparison (|3j;), for 
a chosen Rtf = 100 ^m, which makes it possible to over- 
lay the experimental data (circles) and theoretical curve (line). 
Finally Fig. (Sjl indicates the relative or root-mean-square 
(rms) deviations for these TF fits to theory. This figure was 
made in collaboration with the authors of Ref. This tem- 
perature was chosen to exhibit the TF functional form in the 
regime where it is most problematic, since has a maxi- 
mum there. The experimental data were estimated to corre- 
spond to roughly this same temperature, based on the analysis 
of Ref. 10. Two of the three dimensions of the theoretical 
trap profiles were integrated out to obtain a one-dimensional 
representation of the density distribution along the transverse 
direction: n{x) = J dydzn{r). 

This figure is in contrast to earlier theoretical studies which 
predict a kink at the condensate edg e Moreover, in 

contrast to the predictions of Ref. 1 3, the curves behave mono- 
tonic ally with both temperature and radius. Indeed, in the uni- 
tary regime the generalized TF fitting procedure of Ref. [l^ 
works surprisingly well for our theory with spherical traps, 
and for anisotropic experiments (shown in Fig.|3t). This is, in 
part, a consequence of the fact that trap anisotropy effects be- 
come irrelevant for these one dimensional projections. These 
reasonable TF fits apply to essentially all temperatures investi- 
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FIG. 3: (color online) Temperature dependence of (a) experimen- 
tal one-dimensional spatial profiles (circles) and TF fit (line) from 
Ref. 10, (b) TF fits (line) to theory at T ?a O.TTc ^ 0.1 9Tf (circles) 
and (c) overlay of experimental (circles) and theoretical (line) pro- 
files, as well as (d) relative rms deviations (x^) associated with these 
fits to theory at unitarity. The circles in (b) are shown as the line in 
(c). The profiles have been normalized so that N = J n{x)dx = 1, 



and we set Rtf = 100 fim in order to overlay the two curves, 
reaches a maximum around T — O.lQTp. 
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gated experimentally (loll . and all temperatures we have stud- 
ied, including in the normal state. 

To probe the deviations from a TF functional form, and to 
establish the role of the condensate, in Fig. |5Jl, we show the 
(relative) rms deviation, or from the TF fits as a function 
of T. Except at the very lowest temperatures (where the 
is small and the condensate occupies the entire region of the 
trap), the condensate causes a small, but potentially measur- 
able variation from the TF form, increases rapidly below 
Tc and reaches a maximum around 0.7Tc. Here the profile in- 
volves two different functional forms: the condensate and the 
thermal distribution functions. More precisely, the systematic 
behavior of indicates when T has safely surpassed Tc. 

The reason the Thomas-Fermi fits work well below Tc has 
to do with the presence of noncondensed pairs. If these 
bosonic-like states are ignored, (by omitting the top-most 
shaded regions) in the right panels of Fig. |2j a clear bi- 
modal distribution emerges, as has been predicted in the 
literaturel 1 lllTiil . In this case, TF fits would not work well. 
A residue of this previously discussed bi-modality must be 
present to some extent in the present plots, but this can only 
be seen in the derivatives of our profiles. We find a small kink 
in the first order derivative and a sharp jump in the second or- 
der derivative of n(r), reflecting the condensate edge. These 
features may be more difficult to identify with current experi- 
mental techniques. Alternatively, one may infer the presence 
of a condensate through plots from TF fits. These experi- 
ments will require substantial averaging over a large number 



of profiles, however. 

In summary, in this paper we have found that, in the near 
BEC regime, the component profile for the noncondensed 
bosons, is different from the Gaussian form generally as- 
sumed, although the general shape is consistent with exper- 
iment. This result which follows because these noncondensed 
states are in equilibrium with a condensate, should have im- 
plications for measurements of the condensate fraction. Near 
unitarity, we have found that, our calculated density profiles 
are consistent with experiment and provide strong support for 
using Thomas Fermi fits to the profiles. One can infer that the 
condensate is seen, not in the profile shapes, but, presumably, 
in their temperature calibration. 
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Grant No. DMR-02 13745. 
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